DAHLGREN  DIVISION 
NAVAL  SURFACE  WARFARE  CENTER 

Virginia  22448-5100 


NSWCDD/rR-96/54 


ANALYSIS  OF  THE  PROBABILITY  DENSITY 
FUNCTION  OF  THE  MONOPULSE  RATIO  RADAR 
SIGNAL 


BY  WINSTON  C.  CHOW 
WARFARE  ANALYSIS  DEPARTMENT 


AUGUST  1996 


Approved  for  public  release;  distribution  is  unlimited. 


®irc 


Qtr, 


19961015  035 


TfflS  DOCUMENT  IS  BEST 
QUALITY  AVAILABLE.  THE  COPY 
FURNISHED  TO  DTTC  CONTAINED 
A  SIGNIFICANT  NUMBER  OF 
PAGES  WHICH  DO  NOT 
REPRODUCE  LEGIBLY. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMBNo.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  search  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  or  any  other  aspect  of  this  collection  of  information,  including 
suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Sen/ices.  Directorate  for  information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  99909- 
4302,  and  to  the  Office  of  Management  and  Budget.  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC  20503. 

T.  AGENCY  USE  ONLY  (Leave  blank)  |  2.  REPORT  DATE  I  3.  REPORT  TYPE  AND  DATES  COVERED 


August  1996 


4.  TITLE  AND  SUBTITLE 


Analysis  of  the  Probability  Density  Function  of  the  Monopulse  Ratio  Radar 
Signal 

6.  AUTHOR(s) 

Winston  C.  Chow 

7  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Commander 

Naval  Surface  Warfare  Center 
Dahlgren  Division  (Code  A52) 

17320  Dahlgren  Road 
Dahlgren,  VA  22448-5100 

SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


5.  FUNDING  NUMBERS 


8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 

NSWCDD/TR-96/54 


10.  SPONSORING/MONITORING  AGENCY 
REPORT  NUMBER 


11.  SUPPLEMENTARY  NOTES 


12a.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 


12b.  DISTRIBUTION  CODE 


13.  ABSTRACT  (Maximum  200  words) 

Monopulse  radars  are  often  used  by  the  U.S.  Navy  Department  to  track  airborne  targets  such  as  missiles  attacking  ships.  In  order 
to  effectively  make  use  of  the  target  position  data  provided  by  the  radar,  an  estimate  of  the  radar  output  accuracy  must  be  known. 
This  report  describes  portions  of  research  performed  at  Naval  Surface  Warfare  Center,  Dahlgren  Division  (NSWCDD)  to  develop 
statistical  models  of  monopulse  ratio  behavior  to  improve  accuracy  analysis  of  the  monopulse  radar. 

The  products  of  this  work  include  a  formula  for  the  distribution  of  the  real  valued  part  of  the  radar  signal.  Also  developed  were  two 
tools  to  compare  this  distribution  with  a  normal  distribution.  With  these  formulas,  we  showed  that,  given  certain  assumptions,  this 
distribution  can  be  approximated  by  a  normal  distribution.  Finally,  this  work  showed  the  validity  of  certain  pre-existing  formulas 
developed  at  NSWCDD.  These  results  are  useful  in  developing  algorithms  for  assessing  the  accuracy  of  target  location  estimates  by 
monopulse  radars. 


14.  SUBJECT  TERMS 

monopulse  ratio,  radar  signal,  probability  density,  normal  distribution,  kurtosis, 
Rayleigh  distribution,  random  disturbances,  random  noise,  signal  mean  and  variance 


15.  NUMBER  OF  PAGES 

55 

16.  PRICE  CODE 


17.  SECURITY  CLASSIRCATION 
OF  REPORTS 

UNCLASSIFIED 

NSN  7540-01-280-5500 


18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

UNCLASSIFIED 


19.  SECURITY  CLASSIRCATION 
OF  ABSTRACT 

UNCLASSIFIED 


20.  UMtTATION  OF  ABSTRACT 


Standard  Form  298  (Rev  2-89) 

Prescribed  by  ANSI  std.  239-18 
298-102 


NSWCDD/TR-96/54 


FOREWORD 


Monopulse  radars  are  often  used  by  the  U.S.  Navy  Department  to  track  low  elevation 
airborne  targets  such  as  missiles  attacking  ships.  In  order  to  effectively  make  use  of  the  target 
position  data  provided  by  the  radar,  an  estimate  of  the  radar  output  accuracy  must  be  known. 
This  accuracy  information  is  necessary  in  order  to  judge  how  much  distance  and  at  what  angle 
the  actual  target  can  be  located  away  from  the  radar  estimated  location.  Such  information  is 
necessary  for  defensive  actions  against  the  target.  Also,  radar  accuracy  information  is  needed  in 
order  to  make  efficient  use  of  the  radar  equipment  itself  since  the  amount  of  time  the  radar 
requires  in  tracking  the  specific  position  of  a  target  increases  with  the  desired  accuracy. 

This  report  describes  portions  of  research  performed  at  the  Naval  Surface  Warfare 
Center,  Dahlgren  Division  (NSWCDD)  to  develop  statistical  models  of  monopulse  ratio  behavior 
to  improve  accuracy  analysis  of  the  monopulse  radar.  The  products  of  this  work  are  mainly  new 
formulas  and  confirmations  of  existing  formulas.  These  results  are  useful  in  developing 
algorithms  for  assessing  the  accuracy  of  target  location  estimates  by  monopulse  radars. 
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helping  to  examine  some  of  the  problems  encountered  during  this  research. 
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CHAPTER  1 

INTRODUCTION 


Monopulse  radars  are  used  by  the  U.S.  Navy  Department  for  tracking  low  elevation 
airborne  targets  such  as  missiles  attacking  Naval  ships.  One  such  missile  in  this  category  is  the 
Antiship  Cruise  Missile.  In  order  to  effectively  make  use  of  the  target  position  data  provided  by 
the  radar,  an  estimate  of  the  radar  output  accuracy  must  be  known.  This  accuracy  information  is 
necessary  in  order  to  judge  how  much  distance  and  angle  the  actual  target  can  be  located  away 
from  the  radar  estimated  location.  Such  information  is  necessary  for  defensive  actions  against 
the  target.  Also,  radar  accuracy  information  is  needed  in  order  to  make  efficient  use  of  the  radar 
equipment  itself  since  the  amount  of  time  the  radar  requires  in  tracking  the  specific  position  of  a 
target  increases  with  the  desired  accuracy. 

As  with  any  signal  processing  instrument,  random  disturbances,  commonly  known  as 
random  noise,  are  generated  during  its  use.  These  disturbances  occur  from  both  sources  external 
to  the  radar  as  well  as  from  the  instrument  itself.  The  random  noises  in  the  radar  signal  form 
errors  in  the  radar  reading  that  constitute  a  significant  source  of  inaccuracy  in  target  position 
estimation.  Andyzing  the  nature  of  this  random  noise  is  useful  in  minimizing  the  impact  of  such 
errors  and  making  efficient  use  of  radar  resources. 

Another  reason  of  why  a  better  understanding  of  random  disturbances  is  useful  is  related 
to  how  radar  tracking  is  performed.  Estimating  algorithms  such  as  the  Kalman  filter  are  used  in 
the  tracking  software  in  order  to  estimate  the  position,  velocity,  and/or  acceleration  of  the  target. 
One  of  the  inputs  of  such  estimators  is  the  measurement  error  covariance.  In  the  past,  there  have 
been  no  adequate  approximations  of  this  quantity.  One  other  application  of  the  study  of  random 
disturbances  is  a  better  understanding  of  how  and  (to  what  extent)  the  errors  caused  by  such 
noise  propagate  through  the  radar  system,  causing  inaccuracy  in  the  final  reading. 

Due  to  the  random  nature  of  this  noise,  deterministic  mathematical  methods  do  not  form 
adequate  models  of  radar  signals.  Hence,  radar  signals  needs  to  be  studied  from  a  statistical  and 
probabilistic  point  of  view.  There  has  been  research  work  that  performed  on  the  statistical  nature 
of  monopulse  ratio  radar  signals.  Some  of  this  work  was  documented  in  References  1,  2,  and  3. 
The  additional  new  research  shown  in  this  present  document  supplements  and  extends  these 
investigations.  The  products  of  this  work  are  mainly  in  the  forms  of  new  formulas  as  well  as 
confirmations  of  existing  formulas.  These  results  are  useful  in  developing  algorithms  for 
assessing  the  accuracy  of  target  location  estimates  by  monopulse  radars. 

In  analyzing  the  statistical  nature  of  a  monopulse  radar,  any  information  on  the  nature  of 
the  joint  probability  density  function  of  the  complex  valued  monopulse  ratio  is  important  for 
specifying  the  statistical  parameters.  This  joint  probability  density  function  for  various  cases  has 
been  derived  in  Reference  1,  the  distribution  of  the  real  part  of  the  monopulse  ratio  was  then  still 
unknown.  Yet,  the  probability  density  function  of  the  real  part  of  the  complex  monopulse  ratio 
along  with  any  other  properties  that  can  be  derived  to  further  describe  this  probability  density 
function  is  useful  in  determining  the  errors  of  the  angular  displacement  of  the  target  being 
tracked  by  the  radar.  The  research  described  in  this  report  helped  to  characterize  the  real  part  of 
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the  monopulse  ratio.  In  this  report,  the  probability  density  of  the  real  part  of  the  monopulse  ratio 
for  the  zero  mean  and  partially  correlated  case  is  derived  based  upon  the  joint  density  described 
in  Reference  1.  In  addition,  two  methods  of  analyzing  this  probability  density  function  to  help 
determine  its  closeness  to  a  normal  (Gaussian)  distribution  are  described.  Furthermore,  a 
comparative  study  is  performed  that  compares  some  of  the  results  given  by  Seifer  (References  2 
and  3)  to  the  results  reported  by  Groves  and  Blair  (Reference  1).  In  order  to  enhance  clarity,  an 
explanation  of  the  equivalence  of  a  complex  valued  random  variable  to  a  bivariate  random  vector 
is  given  in  the  last  section  of  this  report.  The  derivation  of  the  variance  of  the  real  part  of  the 
probability  density  function  for  the  zero  mean  and  partially  correlated  case  is  described  in 
Appendix  A.  This  variance  formula  is  needed  in  one  of  the  methods  for  analyzing  the  closeness 
of  Ae  probability  density  to  the  normal  density. 
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CHAPTER  2 

PROBABILITY  DENSITY  OF  THE  REAL  PART  OF  THE  MONOPULSE  RATIO 


bivariate  density  of  the  monopulse  ratio  with  partial  correlations  is  known 
(Reference  1),  the  density  of  the  real  part  of  the  monopulse  ratio  was  not  yet  derived  In  this 
section,  the  probability  density  function  of  the  real  part  of  the  monopulse  ratio  is  derived  as  the 
marginal  density  of  the  joint  probability  density  function  of  the  real  and  imaginary  portions  of 
the  monopulse  ratio.  The  magnitude  of  the  signal  is  assumed  to  exceed  a  given  threshold  level  in 
order  to  be  considered.  This  is  further  explained  in  the  next  paragraph. 

In  a  monopulse  radar,  the  signal  consists  of  four  beams,  which  is  modeled  as  four  random 
vanables,  X,  Y,  U,  and  V.  Each  of  these  random  variables  is  assumed  to  be  normally  distributed. 
Furthermore,  these  four  forms  of  the  signals  are  transformed  into  the  sum  and  difference 
channels.  Let  S=U+iV  be  the  sum  channel  of  the  monopulse  ratio  and  D=X+iY  be  the  difference 

channel.  Let  P+iQ=R=D/S  be  the  monopulse  ratio  and  p  be  the  correlation  between  x  and  u  (the 
real  parts),  which  is  assumed  to  be  the  same  as  that  of  y  and  v  (the  imaginary  parts). 

Furthermore,  let  Cx  and  Ou  represent  the  standard  deviations  of  x  and  u,  respectively,  and 

y=Ox/Ou.  The  threshold  assumption,  mentioned  in  the  previous  paragraph,  is  that  the  magnitude 
of  the  sum  channel  must  exceed  a  given  positive  constant  Rq  in  order  to  be  considered  a  true 
signal.  In  other  words,  the  assumption  states  that  in  order  to  be  recognized  as  a  legitimate  signal, 
the  signal  must  be  such  that  V(u2+v2)  >Ro  >  o. 


Suppose  that  the  monopulse  ratio  is  zero  mean  and  has  the  following  covariance  matrix; 


X  =  a: 


y 

0 

py 


0  py  0 
2 

y  0  py 
0  10 
py  0  1 


(2.1) 


Then,  the  joint  probability  density  function  of  the  real  and  imaginary  parts  of  the  mononulse  ratio 
IS  shown  in  Equation  (7.12)  of  Reference  1  to  be  as  follows: 


r  f  1  y^i-P^h^p2  + 

Pp,Q(P.q)  = - [ - —2  2  2 

"  (p  +q  +y  -2pyp) 


1-11 - g_pyp)  +  lj^-h^(p^+q^+py2-2pyp)  2) 
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where  h  is  a  "normalized"  threshold,  defined  as 


h  = 


«0 

a^^|2(l-p) 


(2.3) 


and  Rq  is  the  actual  threshold  level.  For  the  sake  of  completeness,  one  may  note  that  a  general 
probability  density  function  for  P  and  Q  with  any  given  correlation  of  the  noises  and  without  the 
threshold  restriction  was  derived  and  given  in  Equation  (6.12)  of  Reference  1.  Definitions  of 
various  variables  in  this  equation  are  given  by  Equations  (6.2),  (6.6),  and  (6.13)  of  Reference  1. 
This  formula  for  the  general  probability  density  function  for  P  and  Q  was  later  empirically 
confirmed.  The  empirical  confirmation  was  done  by  programming  the  formula  using  the 
MATLAB  computer  language  and  numerically  integrating  the  formula  over  all  p  and  q.  (The 
numerical  result  of  this  integration  was  one.  See  Appendix  B  for  more  details.) 

The  probability  density  of  the  real  part,  P,  of  the  monopulse  ratio  can  be  derived  as  the 
integral  over  the  real  line  of  the  joint  density  with  respect  to  q.  Hence, 


Pp  (p)  = 


(2.4) 


Y^(l-P^)  y  jh^(P^  +q^  +y^ 


2  2  2  2 
(p  +  q  +  Y  -2pYp) 


2pYp)+l  2.2  2  2 

]exp[-h  (p  +q  +PY  -2pYp)]dq 


Let  ci=p2+72-2p7p  and  c2=p^+p'hF'-2pyp.  Since  by  definition  of  correlation,  the  magnitude  of 
p  is  no  greater  than  one  and  p2+Y2^yp^  one  can  conclude  that  ci=p2+'y2_2pYp^.  Therefore, 

Cl,  being  a  non-negative  number,  can  be  represented  as  ci=c2  for  some  quantity  c.  Using  these 
symbols. 


Pp(p)  = 


Y^(l-P^)  7  h^(q^  +c^)+l  2  2 

J  I - 2 - - ]exp[-h  (q  +  C2  )]dq 

n  ~oo  (q  +  c  ) 


(2.5) 


2  2 

Y  (1-p  )  2  ‘I 

-exp{-h  C2)[  J 


h“  -h^q^  1  _u2„2 


n 


2  2 
-<»  (q  +  c  ) 


dq  +  / 


9  .1  1 

2  2  2  * 


-00  (q  +  c  ) 
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One  can  see  that  this  probability  density  contains  the  sum  of  two  integrals.  Consider  the  first 
integral.  Looking  at  Equation  (3.466)  of  Reference  4, 


where 


2  2 


2  2 

—  cx>  (q  +  c  ) 


-h  q  .  ^.2  f 

dq  =  2h  J 


2  f  e 


1,2  2 

-h  q 


2  2 

2  ^  h  c 

=  2h  [l-(|)(ch)]  —  e 

2c 


0(q^  +  ) 


-dq 


(2.6) 


2  ch  _^2 

<|>(ch)  =  -pr  J  e  dt  => 

Vn  0 


/ 


2  2 


-h2q2  h2^2 

’  dq  = - e 

c 


[I-<|>(ch)] 


— oo  (q  +  c  ) 

Now  consider  the  second  integral  in  the  formula  for  Pp(p).  Note  that 

2 


d  «  h 

^  "2  2 
dc  —oo  q  +  c 


u2^2  « 

e  ^  dq  =  -2h  c  J 


1 


2  2  2 
-oo(q^  +  c  ) 


u2  2 

-h^q^  _ 
e  dq 


(2.7) 


2  2  2 
-cx>(q^  +  c  ) 


—1  d 

e  ^dq  =  (— ^)  — (J 


2  ^  V  2  2 

2h  c  dc  — oo  q  +  c 


1,2  2 

e  ^  dq) 


-1  d  jih^  u2  2 

=  (-^-)  —  l - e** 

2h  c  dc  c 


(l-<t>(ch))] 


(2.8) 


-71  d  e 
=  ( - )  — [- 


1.2^2 

n  c 


2  c  dc  c 


(l-(t)(ch))l 
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Taking  the  derivative  of  the  first  factor. 


2  2 
d  exp (he) 


2  2  2  1 


—  ( - )  =  exp(h  c  )(2h  - -rr) 

dc  c 


2  2 

d  exp(h  c  )  2  2  2  1 

—  [ - (1  -  <|)(ch))]  =  exp(h  c  )(2h  -  — )(1  ~  (t)(ch)) - 

dc  c 


(2.9) 


2  2 

exp(h  c  )  d 


<t>(ch) 


c  dc 

With  the  derivations  just  given,  one  can  conclude  that 

2  2 

7  e  ^  -n  h2  2  2  1  Id 

i  — — ry****^ — ®  -^)(i-(i)(ch))- — (|)(ch)] 

—00  (q+c)  2c  c  ede 


-n  1,2  2  2h^c^-l  d 

—  e  [( - )(l-<l>(ch)) - (t){ch)] 

2c  c  dc 


where 


(2.10) 


d  d  2  cj>  *2  2h  _c2i,2 

—  <|)(ch)  =  — (— )  I  e*  dt  =  — e  ' 
dc  dc  ^ln  0  Vre 


Note  that  the  formula  for  the  derivative  of  (|)(ch)  was  derived  by  using  the  Leibniz  theorem.  Now 
putting  the  two  integrals  together. 


h-  -h2n2  «  1 


2  2  *  dq+  I 


-c»  (q  +  c  ) 


u2  2 

“h  , 

2  rT* 


-00  (q  +  c  ) 


(2.11) 


h^c^  2  2 

Jte  2  1  2h  c  -1  d 

- {h  [1  -  (t)(ch)]  -  —  [( - )(1  -  <|)(ch))-  —  <t)(ch)]} 

c  2c  c  dc 
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Now  recalling  the  equation  for  Pp(p)  given  by  Equation  (2.5),  the  following  must  hold  true: 


Pp{p)  = 


Y  (1-p  ) 


{h  [l-<|>(ch)]- 


1  2h^c^-l  d 

- [( - )(1  -  <|)(ch))-  — (t)(ch)]} 

2c  c  dc 


Y  (1-p  )  hV  1  2h^c^-l 

— - e  e  ^{b  [l-<|>(ch)] - [( - 

c  2c  c 


)(1  -  <t>(ch))] 


(2.12) 


1  2h  _-2,^2 
+  (— )(-p)e  **  } 

2c  yn 


„  ,  ,  Y^(l-p^)  h^c^  -*»%  2  1  2h^c^-l 

Pp(p)  = - e  e  2  {h  [l-(t>(ch)l - [( - 

c  2c  c 


)(l-<t)(ch))]} 


+ 


Y^(l-p^) 


(2.13) 


with  c2=p2+p2^-2pYp  and  c2=ci  and  ci=p2+Y2-2p7p,  as  were  previously  defined.  So  the 
formula  for  the  probability  density  function  of  the  real  portion,  P,  of  the  monopulse  ratio  has  now 
been  derived. 


This  formula  will  now  be  modified  to  a  form  that  can  be  more  easily  programmed  on  a 
computer  for  use  in  monopulse  ratio  radar  signal  processing.  The  given  formula  for  Pp(d) 
implies  that 


Pp(p)  =  Lii_E_i^hV(l  p2)^f^2 


2  2 
2h  c  -1 


2c 


-)][l-(|)(ch)]} 


+ 


Y^(l-P^) 
c^  Vre 


(2.14) 


2-5 


NSWCDD/TR-96/54 


This  implies  that 


Pp(p)  = 


2c^ 


[l-(t)(ch)]e 


hV(l-p^) 


'4n 


he 


■h^c. 


(2.15) 


n  /  N  ^  (1“P  ),  h^Y^(l-D^)  2ch  -h^c, 

•■•Pp(p)  = - ^ - {[1 -(t>(ch)]e  ^  P  ^  +-p-e 

Vrc 


2c 


(2.16) 


Notice  that  the  probability  density  function  of  P  is  obviously  not  normally  distributed.  In 
the  next  two  sections,  formulas  will  be  derived  that  can  be  used  to  provide  some  form  of 
measuring  how  close  the  real  part  of  the  monopulse  ratio,  whose  probability  density  function  was 
just  derived,  is  to  a  Gaussian  density  function.  The  first  method  exploits  the  fact  that  any 
Gaussian  distribution  has  a  kurtosis  value  of  3,  and  the  second  method  uses  a  stochastic  distance 
measure. 
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CHAPTER  3 


FOURTH  MOMENT  ABOUT  THE  MEAN  AND 
KURTOSIS  OF  THE  REAL  PART  OF  THE  MONOPULSE  RATIO 


One  method  of  giving  a  measure  of  the  closeness  between  the  real  part,  P,  of  the 
monpulse  ratio  and  a  normally  distributed  random  variable  uses  the  fact  that  the  kurtosis  of  any 
Gaussian  (that  is,  normal)  distribution  has  a  is  3.  Hence,  finding  the  kurtosis  of  P  and  comparing 
it  with  3  is  one  method  of  comparison  of  the  closeness  of  P  to  a  Gaussian  distribution.  Certainly, 
there  can  be  distributions  other  than  Gaussian  that  also  have  a  kurtosis  of  3,  but  this  is  still  one 
form  of  comparison  between  P  and  a  normally  distributed  random  variable.  Furthermore, 
deriving  the  formulas  for  the  fourth  moment  of  P  along  with  its  kurtosis  may  themselves  be 
useful  information  concerning  the  statistical  nature  of  monopulse  ratio  radars. 

Recall  that  kurtosis  is  defined  as  in  the  following. 


Definition:  Let  X  be  a  random  variable.  The  kurtosis  of  X  is  given  as 


Kurtosis(X)  = 


E[(X-E[X])'*] 

4 


(3.1) 


In  order  to  simplify  notations,  let  t=P-E[P].  Also,  let  Ei  be  defined  as  the  following 
function  for  any  given  real  number  y: 


Ej  (y)  =  I - du  (3.2) 

y  « 

Also,  let  Y  and  p  be  as  previously  defined.  One  may  note  the  fact  that  according  to  Reference  1 
(on  the  top  of  page  7-3), 


E[P]=py  (3.3) 

The  fourth  moment  of  P  about  the  mean,  py,  will  now  be  derived.  With  the  notations  that 

were  just  given,  the  fourth  moment  about  the  mean  can  be  denoted  by  Elt^].  In  finding  the 
fourth  moment,  a  transformation  from  the  Cartesian  coordinates  to  the  polar  coordinates  will  be 
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made.  Hence,  let  t=P-E[P]=Rcos0.  From  Equation  (7.16)  of  Reference  1  the  variance  (second 
moment  about  the  mean)  of  P  is  given  by  the  following  equation: 

E[t^]  = 


(3.4) 


Y^(l-p^)27t  2  7  3  h  [R  +7  (1-p  )]  +  l  2  2 

-  /  cos  Gde/R-^dRC - - - - - - ]exp[-h^R^] 

Jt  0  0  [R^  +Y^(l-p^)]^ 


Writing  the  new  variable  t  in  the  form  of  t=Rcos0,  the  above  formula  implies  that  E[t4]  can  be 
represented  as  follows: 


E[t^]  = 


y^(l-p^)2n  4  oo  2  ^  2  2 

-  i  cos^  0d0/R^dR[ - - - - 1;— - ]exp[-h^R^] 

no  0  [R^  +y^(l_p2)]2 

Looking  at  the  first  integral  in  this  equation,  the  following  relationships  must  be  true. 


(3.5) 


2tc 

/  cos**  0d0  =  [— 0  + 

0  8  4 


4  3  sin 26  sin 46  3 

^  r\  -an  r  ^  .  ..  ^iV 


Iq  =-(27t)  =  -JI 
32  8  4 


(3.6) 


^,4  3  2  2  7  5  h^[R^  +Y^(1-P^)]+1  22 

E[t  ]  =  -Y  (1-p  )/R  dR[ - - - - - - — - - ]exp(-h^R^) 

4  0  [R2+y2^j_p2^j2 


(3.7) 


In  performing  this  integral,  another  change  of  variables  will  now  be  made.  Let  u=r2.  This 
means  that  du=2RdR  and  RdR=du/2.  Substituting  these  values. 


,,,  4  3  2  2  7  2  [u  +  Y^  (1  -  P^  )]  +  1  2 

E[t  ]  =  — Y  (1-p  )Ju  du[ - - - - — - - ]exp(-h  u) 

8  0  [u  +  Y  (1-p  )] 


(3.8) 


3-2 


NSWCDD/TR-96/54 


Now  let  a=u+Y2(l-p2).  This  implies  that  da=du  and  u=a-y2(l-p2).  Substituting  these  values. 


E[t*]  = 


^2  2T  2  22h  a  + 1  2  2  1  2 

y  (1  -  p  )  /[a-y  (1  -  p  )]  da( - r — )exp(-h  a)exp(h^Y^  (1  -  )) 

8  y2(i-p2)  a 


-2  2,  h-y-(l-p-) 


2„2.,  „2.  CO  ^  2,,2  ,  h^a  +  1  -h2a 


=  -y  (l-p  )e 

8  y2(l-p2) 


J[a-y  (l-p^)]^da( - ^ — >« 


a 


2,2  _  p2  ^^h2,2(i_p2)  7j„2  _  ^^^2  _  p2  ^  ^  ^4  _  ^2  ^2  , 

*  y^(i-p^) 


(3.9) 


ha  + 1  _i,2, 


( - 5— )e 


h^a 


a 


In  order  to  further  simplify  formulas,  let 

K=(3/S)'^(  l-p2)exp(h2Y2(i  .p2)) 

Thus, 


(3.10) 


^2~  ^  ^  CO  h^a  +  1  _1,2^ 


Bit'll  =  K{  J(h^a  +  l)e  “da-2Y^(l-p^)  7  (' - )e  “  “"da 

^  ^  y2(i-p2) 


y^a-ph 


(3.11) 


4  22  7  ha  +  1  _h2a 

+y  (I-P  )  J  ( - r — )e  da} 

yhi-ph 


.•.E[t‘‘]  =  K{  Th^oe  ''^“da+  7e  **^“da  -  2y^  (1  -  p^  )h^  7e“‘'^“da 


y^ii-ph  y^dV) 


y^d-P^) 


-2y^(1-P^)  7  (-)e  ‘'^“da  +  y'^  (1  -  p^  7  (— )e“‘*^“da 

y2(l-p2)  a 


2n-n2'»  a 


y^d-p^) 


+y'‘(l-p^)^  7  (^)e  ''^“da} 

y^lV)  “ 


(3.12) 
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f.2  -h'^a,  f  -h'^a.  ^2,,  2..  2  r  -h‘‘a, 

.•.  E[t  ]  =  K{  J  h  ow  da  +  Je  da-2y  (1-p  )h  Je  da 


y^a-ph 


y^a-ph 


y^a-ph 


-»  2,,  2,  7  1,  -h^a.  4,,  2,2  7  h  -h^a^ 

-2y  (1-p  )  J  (— )e  da  +  Y  (1-p  )  1  ( - )e  da 

Y^d-p^)  “  Y^dV)  “ 


4,.  2,2  7  /  ^  ^  , 

+Y  (1-p  )  J  (— :r)e  da} 

r2(i-p2)  o" 


i2  00*2 

i^fi.2  f  -h'^a.  _  ^2,,  2,.  2,  r  -h'^a. 

=  K{h  Ja  e  da  +  [1  -  2y  (1  -  p  )h  ]  Je  da + 

Y^(l-P^)  Y^d-P^) 


(3.12) 


-h^a 


-h^a 


Y^(l-P^)[Y^(1-P^)h^ -2]  J  (- - )da  +  y'*  (1  -  P^  7  )da}> 

.»2/'i  „2-v  ct  .,2/i_„2i 


Y^d-P^)  “ 


Y^d-p^)  “ 


where 


oo  .2 

Ja  e  da  = - 2 — (“**  a-1) 


Y^d-P^) 


Y^d-P^) 


2  2  2  « 
=  [h  Y  (1-p  )+l]- 


-hVci-p^) 


(3.13) 


00-2  P 

,  f  -h  a,  * 
and  J  e  da  =  — 


-h^a 


Y^d-P^) 


-h' 


h2Y2d-p^) 


Y^(l-P^) 


(3.14) 


Now  consider  the  integral  term. 


-h2a 


Y2(i-p2)  a 

-h^a 


-)da 


2  -h^a 


Let  u  =  e  and  dv  =  a  da  =>  du  =  -  h  e  and  v  =  -a 


-1 


(3.15) 
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Substituting  and  then  integrating. 


-h^a 


J  (■ 


■)da  =  -a  e 


-1  -h^a 


-h^a 


-  h' 


J  (- 


-)da 


Y^a-p-^)  y'*(l-p^)  a 


-hVdV) 

1 _ _y,2 

2  2.  ^ 

Y  (1-p  ) 


-h^a 


/  ( 

y2d-p2)  a 


■)da 


Now  let  w=h2a.  Hence,  dw=h2da,  implying  that  da=dw/h2.  Then  note  that 


(3.16) 


-hVdV) 


-h^a 


/  (- 


Y  (1-p  )  y2d_p2)  a 


-)da  = 


^-hVd-p^) 

Y^(l-p^) 


~  e"'^  ^-hVdV)  oo 

.  ^  ( - ^)dw  =  — - - - — -h^  J  ( - )dw 

h^Y^(l-p^)  w/h  y  (l_p  )  j,2y2^j_p2j  w 

(3.17) 

1.2  2r-,  2^ 

-h^Y  d-P  ) 

C  2  7.  7. 

=  — 2 - 2~-^  Ejlh^Y  (1-P  )], 

Y  (1-P  ) 

where  Ei  is  a  previously  defined  function.  This  derivation  shows  that 
oo  ^-hV(l-p^) 

,  ^  -  ( - r">‘*“  =  ~2 - ^-h^Ejlh^Y^d-p^)]  (3.18) 

Y^(l-p^)  «  Y  (1-p  ) 


w  e 

Now  consider  the  integral  term,  J  ( - )da.  Again  letting  w=h2a  giving 

Y^(l-P^)  “ 
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dw=h2da,  implying  that  da=dw/h2.  Using  these  tenns. 


-h-^a 


J  (- 


~  e  “  dw 


y2(i_p2)  a 


■)da=  }  ( - 5-)— r 

h^Y^(l-p^)  w  /  h  h 


-w 


J  ( - )dw  =  Ej  Ih^Y^  (1  -  )] 


hVaV) 


(3.19) 


Y^(l-p^)  a 
Putting  everything  together. 


=0  e~** 

J  (- - )da  =  Ei[h^y^(l-p^)] 

2/1  ^2-) 


(3.20) 


E[t'*]  = 

-h^Y^(l-p^)  -h^Y^(l-p^) 

Kh^  [h^Y^  (1  -  P^  )  +  1]^ - 4 - +  K[1  -  2h^Y^  (1  -  P^  )]- - 5 - 


Ky^  (1  -  p^  )[h^Y^  (1  -  P^  )  -  2]Ei  [h^Y^  (1  -  P^  )]  + 


-h^Y^d-p^) 

Ky^  (1  -  p^  )^  ^ 2 (1  -  P^  )^  Ih^Ej  [h^Y^  (1  -  P^  )1 
Y  (1-p  ) 


(3.21) 


2„..  2  -h^Y^d-p^)  Ke  P  ^  Ke 


-h^Y^dV) 


=  Y  K(l-p  )e 


2y^K(1  -  p^  )e  •>  Y  (1  P  )  _  2y^K(l  -  p^  )Ei  [h^y^  (1  -  P^  )]  + 


4  2  2  ® 

y^Kd-p-*)^- 


-h^Y^d-p^) 

Y^(l-P^) 


(3.22) 
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=  2 


Jj2y2^  2^ 

Ky'*  (1  -  ^ ^  -  2Ky2  (1  -  p2  )Ei  [hS^  (1  -  p^  )] 

y  (i-p  ) 


(3.23) 


^-hVaV) 

=  2K{ - T - -  y^  (1  -  p^  )Ej  [h^y^  (1  -  p^  )]} 


(3.24) 


Now  recall  K  was  previously  defined  to  be  K=(3/8)y2(l-p2)exp(h2^(l-p2)).  Substituting  this 
expression  for  K, 


2K{ 


-hVdV) 

®  2 


-y“(l-p^)Ei[h^y^(l-p^)]}  = 


,  2  2  2  -h^y^(l-p^) 

4 

y^(l-p^)Ei[h^y^(l-p^)]} 


(3.25) 


=  (-)r‘' (I - (1  - p^i 


.•.E[t  ]  = 

-y^  (1  -  ){-\  -  y^  (1  -  p^  )Ei  [h^y^  (1  -  p^  } 

4 


(3.26) 


Hence,  the  equation  for  the  fourth  moment  about  the  mean  of  P  has  finally  been  derived. 

Here  are  some  observations  of  this  fourth  moment  formula.  Notice  that  from  definitions 
already  given, 


2  2 


-u 


Ei[h"y"'(l-p")]=  I  - du,  y  =  -^. 


hV(l-p^)  “ 


and 


2  2 

1 

li^  Ra 


(3.27) 
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The  last  two  of  the  above  three  expressions  implies  that 


1  2y^(1-p^)  2 

2  ~  2 

h  rJ 


(3.28) 


Using  the  formula  for  the  fourth  moment  about  the  mean,  E[t4],  the  kurtosis  of  P  can  now 
be  derived,  based  upon  the  definition  of  kurtosis  and  the  variance  of  P. 

The  variance  of  P  is  given  in  Equation  (7.17)  of  Reference  1.  However,  there  is  a  slight 
error  in  that  formula  since  there  should  be  no  negative  sign  in  the  exponential  expression.  The 
corrected  formula  for  the  variance,  using  the  t  to  represent  P-E[P],  as  before,  is 


E[t2]=l!iL-e^exp(h2Y2(l-p2))Ej[h2y2(l-p2)] 


(3.29) 


See  Appendix  A  for  a  derivation  of  this  corrected  result.  Using  the  definition  for  kurtosis. 


Kurtosis(P)  s 


E[(P-E[P])"] 


a 

1 

~2 


(3  /  4)y^  (1  -  p^  )[-^  -  Y^  (1  -  P^  )Ei  [h^Y^  (1  -  P^  ] 

h 


(1/4)y  (l-p^)^e^ 


(3.30) 


-2h2Y^(l-p2) 


Y^(l-p^)  h^Ej[h^Y^(l-p^)] 


2.,  2, 

Y  (1-p  ) 


h^Y^(l-p^)  2  2  2 

e  ^  ^  P  ^Ei[hV(l-p  )] 


} 


(3.31) 


{■ 


^-2hV(l-p")  h\^l-p2)e-‘*V(lV) 


h^Y^(l-p^)  ‘Ej[h^Y^(l-p^)]  Ejlh^Y^d-p^)! 


3^-hV(lV)  ^-hVdV) 


{- 


2  2  2 
h^Y  (1-P  ) 


h^Y^(l-p^)  'Ej[h^Y^(l-p^)l  Ejlh^Y^d-p^)] 


} 


(3.32) 
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This  implies  that 

Kurtosis(P)  = 

_ 3 _  ^-hVdV)  _  ^2^2  _  p2  _  p2  ^3  33^ 

(1  .  p2  ,,l>V(lV)  ‘  eJ  [b^T*  (1  -  )1  ’ 

This  formula  can  be  used  to  compare  with  the  number,  3,  the  kurtosis  of  any  normal 
distribution.  What  follows  are  some  plots  of  the  kurtosis  of  P  versus  different  values  for  h,r,  and 
g,  respectively.  Notice  how  they  compare  with  the  kurtosis  of  normality,  namely  3.  The  kurtosis 
seems  to  come  close  to  3  for  large  values  of  h  and  g,  while  in  the  case  of  the  r,  correlation 
between  X  and  U,  the  kurtosis  is  close  but  not  equal  to  3  for  small  values  of  r.  The  kurtosis 
grows  large  for  small  values  of  g  and  h,  and  it  also  grows  large  as  r  approaches  1.  This  shows 
that  for  large  values  of  h  and  g  and  for  small  values  of  r,  the  kurtosis  of  the  real  part  of  the 
monopulse  ratio  has  a  kurtosis  that  is  close  to  that  of  a  normal  distribution. 
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HGURE  3-1.  KURTOSIS  OF  P  VS.  THE  y  PARAMETER 
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h  parameter 

HGURE  3-2.  KURTOSIS  OF  P  VS.  THE  h  PARAMETER 
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rho  parameter 

nCURE  3-3.  KURTOSIS  OF  P  VS.  THE  p  PARAMETER 
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CHAPTER  4 

COMPARISON  OF  THE  REAL  PART  OF  THE  MONOPULSE 
RATIO  TO  NORMALITY  USING  A  STOCHASTIC  DISTANCE  MEASURE 


Another  method  of  giving  a  measurement  of  the  closeness  between  the  real  part  of  the 
monopulse  ratio,  P,  and  a  Gaussian  distribution  uses  a  "stochastic  distance  measure"  between  the 
density  of  P  and  the  Gaussian  density. 

Given  two  probability  density  functions,  say,  f  and  g,  the  particular  stochastic  distance 
measure  that  will  be  used  here  is  defined  as  a  quantity,  y,  such  that 


y=  |[f(p)-g(p)]^dp  (4.1) 


For  the  case  being  considered,  let  f(p)=Pp(p),  the  probability  density  function  of  the  real  part  of 
the  monopulse  ratio,  whose  formula  has  already  been  derived  earlier  in  this  paper.  Furthermore, 

let  g(p)=n(p;|i,a2)  the  normal  distribution.  Recall  from  any  statistics  or  probability  theory 
textbook  that 


n(p;p,o  )  = 


l(p-u) 
■exp( - - - ) 


a^|2n  2 

Also,  recall  that  for  the  real  part  of  the  monopulse  ratio,  the  mean,  p,  is  equal  to  py. 


(4.2) 


Using  all  the  formulas  mentioned  above  and  leaving  the  variance,  a^,  as  an  unknown 
quantity,  one  can  evaluate  the  stochastic  distance  measure  integral.  Since  the  stochastic  distance 
measure  must  be  numerically  approximated  on  the  computer,  the  lower  and  upper  limits  of 
integration  needs  to  be  replaced  by  two  large  numbers,  say  L  and  U.  Hence,  the  numerical 
approximation  of  the  stochastic  distance  measure  has  the  following  formula: 

U  2  2 

y(o)=  J  [N(p;pY,a  )  -  P_  (p)]^  dp  ;  L,  U  »  0  (4.3) 

-L 

With  this  construction  in  mind,  a  numerical  minimization  algorithm  can  be  used  to 
approximate  a  a  such  that  y(a)  is  minimized.  This  was  done  by  programming  y(a)  in  the 
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l^^guage,  and  a  MATLAB  routine  known  as  FMIN  was  used  to  estimate  o  such  that 
y(o)  was  minimized.  The  result  was  as  follows: 

Given  inputs:  y  =  l.,  p  =  0.5,  h  =  6.,  L  =  10,  and  U  =  12 

7  (4  4) 

Result:  0  =  0.1157  and y(o)  =  3.6096x10’' 


Thus,  for  Ae  a=0.1 157,  which  was  Ae  approximate  value  Aat  minimized  Ae  stochastic  distance 
measure,  Ae  value  of  Ae  distance  measure  is  as  shown.  The  smallness  of  Ais  number  leads  one 

to  intuitively  believe  that  for  a=0.1157  and  Ae  given  parameters,  PnCp)  can  be  approximated 
accurately  wiA  a  Gaussian  density  function. 
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CHAPTER  5 

SOME  COMMENTS  CONCERNING  FLUCTUATING  TARGETS 


In  the  next  three  sections  to  follow,  a  comparison  is  made  between  some  of  the  results 
found  m  Reference  1  and  those  of  References  2  and  3.  In  Reference  2,  two  cases  of  the  target 
were  considered,  namely  a  nonfluctuating  target  scenario  and  a  target  whose  signal,  other  than 
disturbance  noise,  fluctuates  in  accordance  to  a  Rayleigh  probability  distribution.  The 
comparative  analysis  to  follow  this  section  mainly  deals  with  a  Rayleigh  fluctuating  target.  A 
review  of  the  characteristics  of  the  Rayleigh  distribution  can  be  found  in  Appendix  C. 

According  to  Reference  2  (page  626),  for  a  Rayleigh  fluctuating  target,  E[S]=E[D]=0 
TJis  was  one  of  the  assumptions  for  the  thresholding  case  in  Reference  1.  Also,  as  will  be 
shown  in  the  next  three  sections  to  follow,  the  results  for  the  variance  and  mean  in  the 
thresholding  case  in  Reference  1  are  equivalent  to  those  found  in  References  2  and  3  for  the  real 
part  of  the  monopulse  ratio.  Therefore,  this  implies  that  for  the  thresholding  case  with  the 
covar^ce  matrix  as  previously  given,  the  real  part  of  the  target  signal  other  than  the  noise 
disturbance  of  the  monopulse  radar  in  Reference  1  must  implicitly  also  possess  a  Ravleieh 
distnbution.  ® 

Note  that  in  References  2  and  3,  the  monopulse  ratio  was  defined  only  as  the  real  part  of 
the  complex  valued  monopulse  ratio  as  defined  in  Reference  1.  Also,  this  real  part  is  represented 
by  R  in  References  2  and  3,  while  it  is  represented  by  P  in  both  Reference  1  and  this  present 
report.  Hence,  the  comparative  study  in  the  next  three  sections  to  follow  shall  only  analyze 
properties  of  P,  not  Q.  j  j 
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CHAPTER  6 

COMPARISON  OF  THE  SIGNAL  MEAN  AND  VARIANCE  EQUATIONS 


In  this  section,  the  equations  for  the  mean  and  variance  of  the  real  part  of  monopulse  ratio 
with  thresholding  as  given  in  Reference  1  will  be  compared  to  those  given  in  References  2  and  3. 
These  formulas  will  be  shown  to  be  the  same  as  long  as  any  assumptions  given  in  the  references 
are  properly  taken  into  consideration. 

As  already  stated,  let  D=X+iY  be  the  difference  channel,  S=U+iV  be  the  sum  channel 
and  P+iQ=D/S  be  the  monopulse  ratio.  For  simplicity  in  notation,  let  the  expection  symbol  for 
any  random  variable,  say  A,  be  expressed  as 


E[A1^A  (6.1) 

Finally,  let  the  symbol.  A,  represent  the  event  that  the  signal’s  amplitude  exceeds  the  threshold, 
represented  by  Rq. 

In  Reference  1  (near  the  top  of  page  7-3),  the  mean  of  the  real  part  of  the  monopulse  ratio 
conditioned  upon  the  event.  A,  is  given  by  the  following  equations: 


E[PIA]=: 


®real(D) 

- P  = 

®real(S) 


Ve[(X-X)^1 

Ve[(U-U)^] 


According  to  Reference  2,  Equation  (40), 


E[PIA]  =  E[R]  =  p(a/b) 

where  2a^  sE(ID- Dl^  ]  and  2b^  sE[IS-SI^  ] 
=^E[PIA]  =  E[R]  =  p— ^ 


(6.2) 


(6.3) 
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This  implies  that 


EtPIA]  =  p 


■\/e[IX  +  iY-(X  +  IY)l^  ] 
)/E[IU  +  iV-(U  +  iV)l^  ] 


(6.4) 

_  'y  _  1  1  /'X 

{E[(X  -  X)^  ]  +  E[(Y  -  Y)*^  1}^^^ 

~P  _  2  _  '1  HI  ' 

{E[(U-  U)^]  +  E[(V- 

Keep  in  mind  that  in  deriving  this  equation,  an  assumption  was  made  that  the  amplitude  of  the 
signal  from  the  target  fluctuated  randomly  with  a  Rayleigh  probability  distribution  function. 
Also,  according  to  the  analysis  in  the  previous  section,  the  thresholding  part  of  Reference  1, 
which  is  the  part  of  that  report  being  analyzed  here,  must  also  have  a  Rayleigh  fluctuating  target. 

Note  that  in  Reference  1,  an  assumption  was  made  that  the  variances  of  the  real  and 
imaginary  parts  of  the  sum  channel  are  equal.  Similarly,  the  variances  of  the  real  and  imaginary 
parts  of  the  difference  channels  are  also  assumed  to  be  equal.  In  this  case  the  above  results 
become 


E[PIA]  =  p 


{2E[(X-X)^]}^^^ 

{2E[(U-U)'‘]}^'^ 


(6.5) 


‘^real(D) 

=  P - • 

®real(S) 

This  is  exactly  the  same  result  as  given  in  Reference  1.  Hence,  given  the  assumptions  expressed 
in  the  two  papers  being  compared,  the  equations  for  E[RIRo]  as  given  by  References  1  and  2  are 
identical. 

Now  again  assuming  a  fluctuating  amplitude  target  signal  with  a  threshold  of  Rq,  the 
variances  as  reported  in  the  two  papers  will  be  shown  to  be  the  same.  The  formula  in  Seifer's 
analysis  will  be  given  first  and  then  manipulated  and  compared  with  the  expression  in 
Reference  1.  Assumptions  that  are  used  in  the  derivation  to  follow  will  be  applied. 
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In  Equation  (33b)  of  Reference  2,  the  expression  for  the  signal  variance  conditioned  on 
the  threshold  is  given  by 


2  1  2  *0  'o  (E[IFI])‘^ 

al¥\A]  =  -ii^  exp(-^)Ej  (— ^)  where x  ^  -  -  , 

2  x  +  1  x  +  1  VE[ini  ] 

Ro 

1q  =  normalized  threshold,  -  ~ ,  t]  =  sum  channel  signal  noise, 

VE[ItiI^  ] 

71  _t  (6.6) 

Ej  (x)  s  j-e  dt, 
xt 


a  2  2 

\L  =  —  yl  -  p  ,p  +  i^s  correlation  (D,  S) 

b 

X  as  defined  above  is  interpreted  as  a  "signal-to-noise  ratio." 

Since  a  Rayleigh  distributed  fluctuating  target  is  assumed,  the  sum  channel  signal  can  be 
considered  to  be  the  sum  of  the  noise  plus  the  randomly  fluctuating  actual  sum  channel  signal 
Sinformation  itself,  F.  Hence,  S=F+ti,  which  is  shown  in  Equation  (28)  of  Reference  2.  Note 
that  must  also  be  F+t|  in  Reference  1  for  the  thresholding  case  with  the  given  covariance  matrix 
due  to  the  argument  in  the  previous  section.  According  to  Equation  (31)  of  Reference  2,  the 
following  relationship  can  be  shown: 

2  Eo<Z  +  l)  2  -2 

Iq  = - 5 - ,  recalling  2b  =  Elis -SI  ]  (6.7) 

2b 

2b2  can  be  shown  to  be  equivalent  to  the  sum  of  the  variances  of  the  real  and  imaginary  parts  of 
the  sum  channel  in  the  following  way: 

E[lS-Sl^  ]  =  E[l(u  +  iv)-(u  +  iv)l^  ] 

=  E[l(u  -  u)  +  i(v  -  v)l^  ] 

=  E[(u-u)^  +(v- v)^]^E[IS-SI^ 


1  2  2 

]  =  a„  +cSy 


(6.8) 


•o 


R 


X  +  1 


Ou  +  Oy 


(6.9) 


2  ^ 
o  [PIA]  =  — : 

2b' 


(l-p^ 


)exp(- 


-)Ei(- 


»0 


(6.10) 


Oy  +  Oy 


Ojj  +  Oy 
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2  „2 

setting  -y  =  — y , 
b  <Js 


2 

o^[PIA]  =  -^(l-p^  -C^)exp(- 

2ai 


(6.11) 


-)Ei  (- 


R 


Cu  +  CTy 


0|1  +  Cy 


The  variance  formula,  which  has  a  slight  error,  is  given  in  Equation  (7.17)  of  Reference  1.  The 
corrected  version  of  this  variance  formula  is  derived  in  Appendix  A  and  is  as  follows: 


Var(P)  =  [RIA]  =  (1  -  [h^y^  (1  -  p^  )], 


Ri 


where  h  = 


a/2^ 


P^) 


(6.12) 


The  h  can  be  manipulated  as  follows: 


2  2  2 
h  y^(l-p^)  = 


2  _2  2 
0  2  _  2  _ 

— ^[Y  (1-P  )]  =  7-y-2 


r; 


20x(l-p‘) 


2®x®u 


recalling  that  y  s  , 


12L 


(6.13) 


:.a  [PIA]  =  — y(l-p  )exp( — 5-)Ei  ( — —) 

2^;  2^1  24 

Since  X  and  U  are  the  real  parts  of  the  difference  and  sura  channels  respectively. 


Ox=<Jreal(D)  and  0u=areal(S)  (6. 14) 

Now  a  comparison  between  the  conditional  variance  formula  from  Reference  2  and  the 
modified  formula  for  the  conditional  variance  from  Reference  1  can  be  made.  In  order  to  make 
this  comparison,  one  needs  to  recall  from  a  previous  derivation  that 

=EIIS-SI^  ]  =  a5 +aj  (6.15) 
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By  an  identical  argument. 


Op  =E[ID-DI^  ]  =  ax +Oy  (6.16) 

In  Reference  1,  another  assumption  was  that  the  real  and  imaginary  parts  of  the  variance  of  D  are 
the  same.  Similarly,  the  real  and  imaginary  parts  of  the  variance  of  S  are  assumed  to  be  the  same 
in  this  report.  This  implies  that 

ajy  =2al  and (6.17) 


Furthermore,  Equation  (5.1)  of  Reference  1  indicates  that  correlations  between  parts  of  S  and  D 
do  not  contain  any  imaginary  parts.  Instead,  the  correlation  between  the  real  part  of  S  and  the 
real  part  of  D  as  well  as  the  correlation  between  the  imaginary  part  of  S  and  D  are  represented  by 
the  same  number,  a,  and  cross  correlations  between  real  and  imaginary  parts  are  assumed  to  be 
zero.  These  assumptions  imply  that  which  is  the  imaginary  part  of  correlation  (S,  D),  must 
be  0.  Substituting  all  these  results,  the  two  equations  for  a^^lA]  show  that  the  conditional 
variances  given  in  Reference  2  is  equivalent  to  the  corrected  form  of  the  conditional  variance  in 
Reference  1  provided  that  the  given  assumptions  are  properly  applied. 

This  derivation  is  actually  more  general  than  necessary  since  E[S]  and  E[D]  are  known  to 
be  zero  for  a  Rayleigh  fluctuating  tarjget.  However,  this  result  was  not  substituted,  which  would 
only  slightly  simplify  the  derivation,  in  order  to  avoid  confusion  on  the  meaning  of  the  variances. 

As  a  byproduct  of  this  derivation,  notice  that  given  the  above  stated  assumptions. 


2  2  2 
h^y'‘(l-p^)  = 


X  +  1 


2  2  2 
since  h  y  (1  -  p  )  = 


Ro 


2a 


X  +  1 


(6.18) 


given  that  (y'i)=(Tv  This  result  will  be  used  in  the  analysis  of  the  conditional  variance  with  the 
presence  of  noise  jamming. 
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CHAPTER? 

COMPARISON  OF  THE  SIGNAL  MEAN  AND  VARIANCE 
EQUATIONS  IN  THE  PRESENCE  OF  NOISE  JAMMING 


Another  scenario  considered  in  Reference  2  was  the  case  where  noise  jamming  also 
disturbed  the  radar  target  signal.  This  noise  represents  an  intentional  attempt  to  deceive  the  radar 
from  accurately  tracking  the  target. 

Several  assumptions  were  made  for  the  analysis  in  reference  to  the  jamming  by  Seifer  in 
Reference  2  (page  633).  These  assumptions  are  as  follows: 

•  jamming  comes  from  a  single  point  source. 

•  the  single  point  source  is  within  the  principal  lobe  of  the  receiving  antenna's 
sum  pattern  and  also  copolarized  with  the  antenna  response  there. 

•  jamming  dominates  thermal  noise  of  the  receiver. 

The  analysis  of  the  means  of  the  two  papers  will  now  be  given.  As  previously  given,  the 
formula  from  the  conditional  mean  for  a  Rayleigh  fluctuating  target  according  to  Equation  (40) 
of  Reference  2  is 


EIPIA]=E[P]=p-  (7.1) 

b 

while  the  same  expectation  according  to  Reference  1  (near  the  top  of  page  7-3)  is 

E[piA]=^ri5KDip  (7.2) 

°real(S) 

such  that,  as  before,  p  is  the  real  part  of  correlation  of  D  and  S  for  both  papers. 

Let  Dj  be  the  complex-valued  quantity  representing  the  contribution  on  the  difference 
channel,  and  let  Gj  be  the  same  on  the  sum  channel.  Let  rj  be  Dj/Gj.  According  to 
Equation  (72)  of  Reference  2, 
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E[PIA]  =  p-  = 
b 


X  +  1 


(7.3) 

E[PIA]-f  =  - 

X  +  1 


where  f  is  defined  as  the  real  part  of  the  monopulse  ratio  without  random  disturbance.  Since  the 

meaning  of  a  and  b  is  the  same  as  previously  defined,  and  the  equations  for  E[RIA]  in  the  two 
reports  are  equivalent  under  the  assumptions  given,  this  equation  for  the  jamming  case  must  also 
apply  to  the  model  as  given  in  Reference  1. 

The  variance  equations  for  the  two  reports  under  noise  jamming  will  now  be  analyzed. 
Recall  that,  with  the  given  assumptions, 

2  /  /I  2  2 

H  -C  j  > 


+  a. 


2a 


+  a. 


(1  -  p"-  )  =  —f(l  -  P^  )  =  (1  -  p^  ) 

2o„ 


(7.4) 


1  2 
2 


1  2  2 
=  -Y^(l-p^) 
2 


The  result  as  shown  in  Equation  (72)  of  Reference  2  implies  that 


1 

2 


2 


1  X 

2  (X  +  l)^ 


(7.5) 


1  2  2  ^  1 
.-. -Y  (1-p  )  =  - 
2  2(x  +  l) 


X  ,  .,2 

— '••j  -rl 


Substituting  this  result  into  the  equation  for  the  conditional  variance  according  to  Reference  1 
gives 


a^[PlA]  = 

fX  ^  2  2  2  2  2  2  2 

- ^Irj-rl  exp(h  Y  d  -  P  )Ei  (h  y  (1-p  )) 

2  (X  +  1) 
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(7.7) 


which  is  the  same  as  the  formula  in  Equation  (76)  of  Reference  2. 

Thus,  the  conclusion  of  this  derivation  is  the  fact  that  the  conditional  variances  equations 
of  Reference  1  and  Reference  2  are  equivalent  under  the  given  hypotheses  even  in  the  presence 
of  a  single  point  noise  jamming. 
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CHAPTER  8 

COMPARISON  OF  THE  SIGNAL  MEAN  AND 
VARIANCE  EQUATIONS  FOR  TRACKING  IN  NOISE 


In  July,  1994,  Arnold  Seifer's  correspondence,  supplementing  his  work  shown  in 
Reference  2  appeared  in  Reference  3.  Among  other  new  assumptions,  two  of  them  are  the 
hypotheses  that  no  noise  jamming  is  present  and  random  disturbance  noise  between  the  sum 
channel  and  difference  channel  are  uncorrelated. 

This  lack  of  correlation  between  the  random  disturbance  of  the  sum  and  difference 
channel  signal  does  not  necessarily  mean  that  S  and  D  are  uncorrelated,  but  only  that  the  random 

disturbance  noise  portions  of  S  and  D  are  uncorrelated.  Hence,  p  is  not  necessarily  0  although, 
recalling  one  of  the  hypothesis  of  Reference  1,  the  imaginary  part  of  the  correlation  between  D 

and  S  (namely,  Q  is  assumed  to  be  0. 

The  derivations  in  Reference  3  for  the  conditional  mean  and  the  conditional  variance  are 
based  upon  the  assumption  of  the  lack  of  correlation  between  the  random  disturbance  of  the  sum 
and  difference  channels.  This  assumption  was  used  to  obtain  the  following  relation,  which  can 
be  found  in  Equation  (16)  of  Reference  3 


a  X  ,  « 

—  (p  +  iQ  = - f ;  f  real  by  assumption  (8. 1) 

b  Z  +  1 

Also,  as  Seifer  stated  on  the  same  page,  the  monopulse  ratio  without  the  disturbance  noise  is  by 

definition  real.  Therefore,  this  fact  along  with  the  assumption  in  Reference  1  that  ^=0,  each 
implies  that 

a  X 

-p  = - f  (8.2) 

b  x  +  l 

which  was  also  shown  in  Reference  3  on  the  same  page.  These  two  relations  were  the  only  ones 
that  were  used  in  Reference  3  for  deriving  new  formulas  for  the  conditional  mean  and  variances 

for  the  Rayleigh  fluctuating  case  conditioned  upon  the  event,  Rq. 

Since  the  assumption  of  0  correlation  between  the  random  disturbances  of  D  and  S,  along 
with  the  above  equations  derived  from  this  assumption,  does  not  violate  any  relations  given  in 
Reference  1,  these  results  in  Reference  3  for  the  conditional  mean  and  conditional  variance  of  a 
Rayleigh  fluctuating  target  are  still  equivalent  to  the  formulas  derived  in  Reference  1.  As  noted 
above,  jamming  was  not  considered  in  Reference  3. 
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CHAPTER  9 


EQUIVALENCE  OF  A  COMPLEX  VALUED  RANDOM 
VARIABLE  TO  A  TWO-DIMENSIONAL  RANDOM  VECTOR 


This  paper  has  dealt  mainly  with  the  real  part  of  the  monopulse  ratio  as  a  random 
variable.  Combining  the  random  variable,  P,  with  the  imaginary  part,  Q,  the  monopulse  ratio 
becomes  a  complex  valued  random  variable  P+iQ.  Another  obvious  way  of  representing  this  is 
by  using  the  ordered  pair  (P,Q)  or  its  corresponding  random  vector  (the  transpose  of  (P,Q)). 
Although,  intuitively,  complex  and  vector  representations  seem  clearly  statistically  equivalent, 
complex  numbers  are  mathematically  treated  differently  from  ordered  pairs.  This  difference  may 
cause  some  people  to  still  question  whether  the  two  representations  are  equivalent.  Since  the 
difference  in  the  mathematical  treatment  between  complex  numbers  and  ordered  pair 
representations  are  not  considered  in  determining  whether  the  two  entities  are  probabilistically 
the  same,  they  are  in  fact  equivalent.  The  explanation  to  follow  gives  a  justification  for  this 
equivalence.  This  explanation  will  be  in  the  form  of  a  theorem. 

Theorem  1;  Given  a  complex  valued  random  variable,  X+iY,  there  exists  a  two-dimensional 
random  vector,  namely  (X,  Y)^,  (or  bivariate  distributed  ordered  pair,  namely  (X,Y))  that  is 
statistically  equivalent  to  the  complex  valued  random  variable. 

A  proof  of  this  theorem  is  given  in  Appendix  D.  With  this  theorem  in  mind,  we  can  use  the  three 
forms  of  the  bivariate  random  variable  in  X  and  Y  as  given  in  the  statement  of  the  theorem, 
interchangeably. 

For  a  precise  measure-theoretic  definition  of  a  random  variable,  one  is  referred  to  any 
textbook  on  fundamental  probability  theory  that  is  based  upon  a  measure-theoretic  approach. 

Due  to  this  equivalence  that  is  now  justified,  the  random  variable  P+iQ  and  the  ordered 
pair  (P,  Q)  as  well  as  (P,  Q)T  can  be  discussed  interchangeably  as  a  model  of  the  real  and 
imaginary  portions  of  the  radar  monopulse  ratio.  This  should  help  clarify  any  ambiguity. 
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CHAPTER  10 
SUMMARY 


Several  results  have  been  produced  by  the  investigative  project  described  here.  A  closed 
form  formula  of  the  probability  density  function  of,  P,  the  real  part  of  the  monopulse  ratio  of  a 
monopulse  radar  was  derived.  Furthermore,  two  methods  of  comparing  this  probability  density 
function  with  a  Gaussian  density  were  derived.  One  was  based  upon  kurtosis,  and  the  other  one 
was  based  upon  a  stochastic  distance  measure.  A  third  method  that  could  be  used  to  check 
closeness  to  normality  is  to  compare  the  areas  under  the  curve  of  the  tail  regions  of  the 
probability  density  function  of  P  and  a  normal  probability  density  function.  This  can  be 
investigated  in  the  future.  We  showed  that  the  kurtosis  came  close  to  that  of  a  normal 
distribution  for  large  values  of  the  normalized  threshold,  h,  large  values  of  y=<Tx/Ou>  a^^d  small 
values  of  the  correlation  between  X  and  U  (namely  p).  We  also  showed  that  under  certain 
assumptions,  the  stochastic  distance  measure  of  the  normal  distribution  and  the  distribution  of 
the  real  part  of  the  monopulse  ratio  is  small.  A  comparative  analysis  was  performed  in  relation  to 
some  of  the  formulas  for  the  mean  and  variance  of  P  as  presented  in  Reference  1  with  those  of 
References  2  and  3.  We  showed  that  given  the  correct  hypothesis,  the  formulas  were  equivalent 
in  spite  of  the  fact  that  the  equations  had  different  forms  and  were  derived  from  different 
approaches. 

In  conclusion,  a  closed  form  formula  for  the  distribution  of  P,  is  now  derived  and 
available  for  use  in  evaluating  the  accuracy  of  monopulse  radars.  We  also  have  derived  two  tools 
in  comparing  the  real  part  of  the  monopulse  ratio  distribution  with  the  normal  distribution.  Also, 
using  these  tools,  we  showed  that  given  certain  assumptions,  P  can  be  approximated  as  a 
normally  distributed  random  variable.  Finally,  we  showed  that  despite  the  different  approaches 
of  Reference  1  as  compared  to  References  2  and  3,  the  formulas  for  the  variance  and  mean  of  P 
with  and  without  jamming  are  equivalent. 
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APPENDIX  A 

DERIVATION  OF  THE  VARIANCE  OF  THE  REAL 
PART  OF  THE  MONOPULSE  RATIO  WITH  THRESHOLDING 


A-1/-2 


NSWCDDArR-96/54 


In  the  section  on  kurtosis,  the  variance  of  the  real  part  of  the  monopulse  ratio  was  given 
in  Equation  (3.29)  as 


E[t^  ]  =  ^ ^exp(h^Y^  (1  -  ))Ej  [h^y^  (1  -  )] 


(A.1) 


Recall  that  t  was  defined  to  be  P-E[P],  and  E[P]=pY,  Also  recall  that  this  formula  differed 
slightly  from  the  one  shown  in  Reference  1  of  the  main  text.*  What  follows  is  a  derivation  of  the 
formula  for  the  variance  as  given  here.  Recall  from  the  section  on  kurtosis  that  according  to 
Equation  (7.16)  of  Reference  1,  the  variance  can  be  obtained  from  the  following  integral. 


2  2 
2  y  (1  -  p  ) 
E[t  ]  = - - — * 


2jc 


2  7  3  +y^(l-p^)]  +  l  2  2 

J  cos  ede/R  dR[ - z - z - z—z - ]exp[-h  R  ] 

0  0  [R^  +y^(l-p^)]^ 


(A.2) 


2  2  7  3  h^[R^  +y^(l-p^)]  +  l  22 

=  y  (1-p  )jR^dR[ - - - - - r— ; ]exp[-h^R^] 


0 


2  2  2  2 
[R'‘  +y^(l-p^)]^ 


In  performing  this  integration,  we  need  to  do  a  change  of  variable  by  letting  u=r2.  Then 
du=2RdR,  and  RdR=(l/2)du.  Substituting, 


2  1  2  2  7  h^[u  +  y^(l-p^)]  +  l  2 

E[t  ]  =  — y  (1-p  )Judu[ - - - - — - - ]exp[-h  u] 

2  0  [u  +  y2(i_p2)]2 


(A.3) 


Now  changing  variables  again,  let  a=u+g2(l-r2),  implying  that  du=da  and  u=a-g2(l-r2).  This 
gives 


E[t^]  = 


1  2  2  h  T  2  2  h  oc  +  1  — 

-y  (1-p  )e  J[a-y'®(l-p^)]da( - ; — )e  **  “  = 


Y^dV) 


(A.4) 


a 


7  +  K  -h^a.  2.^  2.  7  +  -h^a,  . 

K{  J  ( - )e  da-Y  (1-p  )  J  ( - 5 — 

y2(l-p2)  “  y^(l-p^)  tt 

for  K=(l/2)y^(l-p2)exp(li2y2(i.p2)). 


(A.5) 


*  Groves,  G.  W.  and  Blair,  W.  D.,  Statistical  Studies  of  the  Monopulse  Ratio, 
NSWCDD/TR-94/97,  Oct  1994,  Naval  Surface  Warfare  Center,  Dahlgren  Division, 
Dahlgren,  VA. 
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Multiplying  out  the  terms, 


E[t^]=K{h^  Te  '*^“da+  T 


y2(1-p2) 


y2(i_p2)  a 


Y^d-p^lh^  7 

Y^dV)  “ 


(A.6) 


-Y^(l-p^)  J  (-Y)e  “da} 
Y^dV) 


E[t^]  = 


7ey“da  +  [l-Y^l-p^h^]  7 


Y^d-P^) 

-Y^(l-p^)  7  (-^)e~*'^“da} 

Y^dV)  “ 


Y^d-P^)  “ 


Now  consider  each  integral  term  separately. 


1  2  A 

f  “•*  “J  ® 

J  e  da  =  — 


-h^a 


Y^d-P^) 


h2Y2(i_p2) 


Y^dV) 


Now  consider  the  integral, 


(A.7) 


(A.8) 


1  1-2 
/  (~T^* 

Y^dV)  « 


Integration  by  parts  will  be  used  for  this  integration. 

^  2 
Let w  =  e  and dv  =  a  da.  .'.  dw  =  -h^e  ^  “da and  v  =  -a  ^ 


1  ,  -h^a 


J  “da  =  -a-"e 

Y^d-P^)  “ 


-1  -h-^a 


Y^d-P^)  Y^dV)  “ 


.2  7  l,-h2a^ 

-h  J  (— )e  da 


-h^Y^d-p^)  - 

e  '  '  f'  ^  -  oo  I  2- 

2  2  "**  ^ 

Y  (1-p  )  Y^d-P^)  “ 


(A.9) 


(A.10) 
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Now  suppose  s=h2a.  Then  ds=h2da  and  da=ds/h2.  Substituting, 


-h2y2(i_p2) 


2  ^  ^ 1 ^  ® 


-hVaV) 


2/,  2-  ■»  -> 

Y  (1-p  )  y2(j_p2)  a 


/  (-)e-"“da  = 


Y^(l-P^) 


7  e-*  .  eo 

■y  J  - T)ds  =  -1 - 2— -•> 

h^y^(l-p^)  s/h  y  (l_p  )  h^y^(l-p^)  s 


/  { - )ds 


(A.11) 


-h2y2(i-p2) 

=  — - —  -  h^Ei  (1  -  p^  )] 

Y  (1-p  ) 


1  -h2 


2  2  ^ 
y2(i-p2)  “ 


-hVd-p^) 

h  ot--  *  .  2„  ,.2  2..  2 


da  = 


- 2 - •»  Ejih  y^'d-p")] 

Y  (1-P  ) 


(A.12) 


Finally,  consider  the  integral  term. 


1  t.2 

f  i.-"  “da 

y2(i_p2) a 


Let  t— h2a.  Then  dt=h2da  and  d(x=dt/h2.  This  implies  that 


(A.13) 


1  l2  oc 

f  ^  -h'^a .  f 

)  —e  da=  J 


1  _d  dt 


y2(i_p2) a 


hV(l-p^)  1  / 


(A.14) 
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Inserting  the  formulas  for  these  three  integrals  into  the  equation  for  E[t2], 


2  2  « 
EIt^]  =  K{h^(- 


-h2y2(i_p2) 

Y^d-P^) 


(A.15) 


h^Ej  (h^Y^  (1  -  P^  ))!}  =  K{Ei  [h^Y^  (1  -  P^  )!} 

=  -Y^  (1  -  P^  [h^Y^  (1  -  P^  )] 

2 

Var(P)  =  E[t^  ]  =  iy2  d  _  p2  [h^y^  d  -  P^  )] 


(A.  16) 


This  is  the  closed  form  formula  for  the  variance  of  P.  Keep  in  mind  that  this  is  a  conditional 
variance,  conditioned  upon  the  event  that  the  sum  channel  signal  amplitude  exceeds  the 
threshold. 
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APPENDIX  B 

A  GENERAL  JOINT  DISTRIBUTION  OF  THE  MONOPULSE  RATIO 
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This  document  mainly  dealt  with  the  radar  monopulse  ratio  with  a  specific  type  of 
covariance  matrix  between  S  and  D  and  with  the  condition  of  thresholding.  A  generalized 
bivariate  distribution  of  P  and  Q  was  given  in  Reference  1  of  the  main  text*  (page  6-3).  The 
equations  are 


Ke-^ 

Pp,Q(P.q)  = - 2  I - 

InA  ^J-\I,\A 


,  where 


A  s  -  4AC;  M  s  F  +  A‘^  (AE^  -  BDE  +  CD^  );  K  s  B^D^  +  B^E^ 
-4BCDE  +  4C^D^  -  4ABDE  +  4A^E^  +  8A^C  -  2B^C  +  8AC^  -  2AB^ 

A  _  ^  2  1  2  1 

A=  ajjp  +  aj2Pq  +  aj3p  +  — a22q  +a23q‘‘‘~“33 
2  2  2 

2  2 

B  = -aijpq  +  aj2  (P  -q  )  ~  “13^  +  ai4P  +  a22Pq  *23P  ^24^  *34 

2  2  1  2  1 

C  =  -ajjq  -  aj2Pq  -  ai4q  + -a22P  +a24P+~*44 

2  2  2 

D  s  -a  jjxp  -  aj2  (yp  +  xq)  -  a  J3  (up  +  x)  -  a  j4vp  -  a22yq  -  a23  (uq  +  y) 
-a24^q-*33“-*34^ 


E  =  ajjxq  -  aj2  (xp  -  yq)  +  aj3uq  -  a|4  (x  -  vq)  -  a22yp  -  a23up 
-a24(vp+  y)-  a34U  -  a44V 


F 


1 

2 


-2  _  _  _  1  _2 
11*  +  *l2*y  +  *13*“  +  *14*'' +  “*22y 

2 


+  a23yu  +  a24yv 


+  -*33“  +* 

34  “V  + 

-*44V 

A 

ai2 

®13 

ai4> 

Z”'  = 

ai2 

a  22 

®23 

^24 

®13 

^23 

®33 

^34 

^^14 

^24 

^34 

®44  j 

(B.l) 


*  Groves,  G.  W.  and  Blair,  W.  D.,  Statistical  Studies  of  the  Monopulse  Ratio, 
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This  generalized  bivariate  probability  density  function  was  empirically  tested  by 
numerically  doubly  integrating  the  function  with  respect  to  p  and  q  over  the  Euclidean  plane  to 
see  whether  the  integration  resulted  in  a  value  of  1.  The  input  parameters  were  chosen  to  be  as 
follows: 


s 


-1 


100 

2 

1,5 

3 

2 

100 

-3 

3.5 

1,5 

-3 

200 

4.5 

3 

3.5 

4.5 

200 

X  =  10,  y  =  5,  u  =  12,  V  =  -6 


(B.2) 


The  programming  for  the  numerical  integration  was  performed  in  the  MATLAB  computer 
package  language.  This  numerical  double  integration  indeed  produced  a  result  of  1.  In  other 
words. 


/  J  - 2  "/  dpdq=»l  (B.3) 

-oo-co2nA 

Since  the  input  parameters  were  arbitrarily  chosen,  this  further  supports  the  correctness  of  the 
generalized  probability  density  function  as  derived  in  Reference  1  of  the  main  text. 
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APPENDIX  C 

THE  RAYLEIGH  DISTRIBUTION 


C-l/C-2 
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A  Rayleigh  probability  function  with  a  parameter  of  a  on  a  random  variable  z  is  defined 
on  the  non-negative  real  line  as  follows: 


(z;a)  = 


-^re  Vze[0,oo) 

1  o 


(C.l) 


0  Vz  elsewhere 


This  distribution  is  very  useful  in  modeling  physical  phenomena  since  it  represents  the 
magnitude  of  the  vector  sum  of  two  independent  and  zero  mean  normally  distributed  random 
variables  with  the  same  variance.  This  fact  can  be  more  clearly  stated  when  written  in  a  form  of 
a  mathematical  theorem. 


Theorem  lA:  Let  X  and  Y  be  independent,  zero  mean,  normally  distributed  random  variables 
with  the  same  standard  deviation,  c.  Then  the  transformation 


z  = 


+  Y' 


(C.2) 


has  a  Rayleigh  distribution  with  parameter  or. 


Proof:  Since  X  and  Y  are  independent  normally  distributed  random  variables  with  a  common 
variance  a, 


Px,Y(*>y)  = 


■  20^ 


;  z 


=  Vx^  +Y^  <z=> 


2  TVS 


(C.3) 


.  [2  2 


F2(z)  =  Prob{(x,y);Vx"  +  y"  <z};  s  cumulative 
distribution  function,  Prob  =  probability 


Let  X  =  rcos9,  y  =  rsinS 

1 


2 

,r  =  Vx  +y 


=»Fz(z)  = 


2a^ 


dxdy  = 


2Ka^  0  0 


ZE?. 
z2it  2 

J  I  e  ®  rdrdO 


(C.4) 


(C.5) 


Performing  the  double  integration,  the  cumulative  distribution  function  for  the  random  variable  Z 
is 


Fz(z)  =  -e^<^^  +l=>Pj^(z)  =  -i[F^(2)]^_Le2a2 

dz  o 

This  is  the  formula  for  the  Rayleigh  probability  density  function,  and  therefore,  the  proof  is 
complete. 
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APPENDIX  D 

PROOF  OF  THE  EQUIVALENCE  OF  A  COMPLEX  VALUED 
RANDOM  VARIABLE  TO  A  TWO-DIMENSIONAL  RANDOM  VECTOR 
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The  following  is  the  proof  of  theorem  1,  given  in  this  report.  Recall  that  this  is  the 
theorem  to  claim  that  a  complex  valued  random  variable  is  statistically  equivalent  to  its  random 
vector  and  bivariate  ordered  pair  counterparts.  The  theorem  is  stated  again  here. 

Theorem  1  (Repeat):  Given  a  complex  valued  random  variable,  X+iY,  there  exists  a  two 
dimensional  random  vector,  namely  (X,  Y)T,  (or  bivariate  distributed  ordered  pair,  namely 
(X,Y))  that  is  statistically  equivalent  to  the  complex  valued  random  variable. 

Proof:  Suppose  X  and  Y  are  random  variables.  Let  V  be  a  random  vector  such  that  vT=(X,Y), 
where  the  superscript  T  refers  to  the  transpose  of  a  vector  (matrix).  Then  by  the  mathematical 
definition  of  a  random  vector,  V  is  a  function  that  can  be  written  in  the  following  form: 


V:n->R^ 


such  that  V(a))  = 


S(ffi) 

Q((o) 


V  (0  e  n,  a  sample  space 


(D.l) 


Notice  that  there  exists  a  one  to  one  and  onto  mapping,  namely  g(x,y)=x+iy,  that  maps 
from  the  R^  space  to  the  complex  plane,  C.  Combining  this  map  with  the  random  vector 

mapping  gives  the  composition  mapping,  Z  =  goV:fl  ^  R  ->  C.  With  this  construction 
along  with  the  mathematical  definition  of  probability  on  a  random  variable,  notice  that 

P[x  +  iy  6  A  c  C]  =  P[Z  6  A  C  C]  =  P{co  e  Q:  Z(a>)  e  A  c  C} 

=  P{(0  €  Q:  V(co)  €  g“'  (A)  c  R^  }  =  P[(x,y)  e  g"^  (A)  c  R^  ]  (D.2) 

=  P{(x,y):g(x,y)  =  X  +  iy} 

Pictorially,  this  can  be  seen  in  the  following  way: 


g 


-1 

g 


q 


HGURE  D-1.  ONE  TO  ONE  ONTO  TRANSFORMATION  OF  EVENT  A 


These  arguments  shows  that  probabilistically,  Z=X+iY  and  V=(X,Y)  must  be  equivalent,  or 
analogously,  Z=X+iY  and  the  vector  V"^=(X,Y)"r  must  be  equivalent.  Therefore,  the  ordered 
pair  notation  and  complex  notation  of  a  pair  of  random  variables  are  statistically  equivalent.  This 
completes  the  proof  of  the  theorem. 
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